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Abstract. Multilevel observations 
of formaldehyde (H2CO) and methanol (CH3OH) toward 
the nearby spiral galaxies NGC 253, Maffei2, IC342, M82 
and NGC 6946 are presented. H2CO was detected in all 
galaxies (tentatively in NGC 6946). CH3OH was detected 
in all objects with the notable exception of M82. 

H 2 CO line intensity ratios point out differences in gas 
density both between galaxies and within the central re- 
gions of individual objects. Model calculations show that 
the bulk of the gas emitting H2CO in NGC 253 is at a 
density of ~ 10 4 cm~ 3 , while the H 2 CO lines in M82 and 
IC342 trace two different, spatially separated gas compo- 
nents with densities of < 10 4 cm -3 and ~ 10 6 cm -3 . The 
south-western molecular hotspot in M82 and the center 
of IC 342 are the regions with the highest density. 

Methanol is subthermally excited in all galaxies, with 
the lowest excitation temperatures found in IC342. The 
CH3OH abundance in NGC 253 and the non-starburst nu- 
clei of IC 342 and Maffei 2 are comparable. A map of the 
3k — 2k lines in NGC 253 shows that CH3OH traces clumpy 
structures better than other molecules requiring high gas 
density to be excited. CH3OH toward M82 is at least 
an order of magnitude less abundant than in otherwise 
comparable galaxies. This confirms the existence of global 
chemical differences, and thus very large scale variations in 
the state of the molecular gas phase, even between galaxies 
commonly classified as starburst nuclei. 
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1. Introduction 

The central regions of galaxies, including our own, are 
characterized by a unique population of molecular clouds, 
both denser and warmer than typical disk clouds. In galax- 
ies undergoing a nuclear starburst, these clouds are sub- 
ject to turbulence and intense UV radiation, as well as 
shocks, due to both a steep gravitational potential and the 
formation and death of young massive stars. They form a 
distinct type of interstellar molecular environment, differ- 
ent from Galactic dark clouds and the surroundings of H II 
regions. Galactic center clouds are a rich source of molecu- 
lar emission; ~ 25 molecular species have been discovered 
in the nuclear regions of external galaxies. 

Formaldehyde and methanol are both useful tracers of 
the dense interstellar gas phase. Their molecular struc- 
ture is more complicated than that of HCN or CS, the 
"classical" tracers of dense gas. This leads to a greater 
wealth of transitions which can be observed in the mil- 
limeter wavelength range, providing a powerful tool to 
understand the physical and chemical conditions of the 
gas associated with starbursts. 

The first formaldehyde transitions detected in external 
galaxies were A"-doublet ortho transitions in the centime- 
ter wavelength range (lio — In at 6cm and 2n — 2i 2 at 
2 cm; Gardner & Whiteoak 1974, 1976, 1979, Whiteoak & 
Gardner 1976, Cohen et al. 1979, Graham et al. 1978, 
Seaquist & Bell 1990, Baan & Goss 1992). Galaxies 
where these lines have been detected include NGC 253, 
NGC 4945, M82, CentaurusA, NGC 3628, M31 and the 
LMC. Recently, the 2n — 2i2 transition has been de- 
tected toward the Einstein Ring gravitational lens system 
B0218+357 at a redshift of z = 0.68 (Menten & Reid 
1996). In all the above sources the lines are observed in 
absorption; with T ex < 2.7 K, they can absorb even the mi- 
crowave background. The lio — In line emission was first 
observed toward Arp220 (Baan et al. 1986). Detections in 
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nine other galaxies and an interferometric map of Arp 220 
(Baan et al. 1993, Baan & Haschick 1995) followed. These 
lines are interpreted as H2CO megamasers. 

There are less observations of the rotational transitions 
at millimeter wavelengths, even though these transitions 
are sensitive tracers of cloud temperature and density. The 
3o3 — 2 2 para transition was discovered in M 82 (Baan et 
al. 1990), the 32i — 220 para line was observed in NGC 253 
(Petuchowski & Bennett 1992), and the 3i2 — 2n ortho 
transition was detected in the LMC by Johansson et al. 
(1994) and tentatively in the southern galaxy NGC 4945 
(Mauersberger et al. 1996b). 

The first successful observation of extragalactic 
methanol (CH3OH) was reported in 1987 by Henkel et 
al.. They detected a superposition of four 2k — lk tran- 
sitions at 96 GHz in NGC 253, IC342 and (tentatively) 
NGC 6946. This group of transitions requires only low ex- 
citation and all four components are among the strongest 
methanol lines observed in Galactic sources. The same 
group of transitions was detected in the southern starburst 
galaxy NGC 4945 (Henkel et al. 1990). The 0q-1_i.E tran- 
sition was seen toward NGC 253 by Henkel et al. (1993). 
Methanol masers, very common in the Galaxy, are also 
detected toward the Large Magellanic Cloud at 6.6 GHz 
(Sinclair et al. 1992, Ellingsen et al. 1994). Here, we re- 
port observations of six H2CO and seven CH3OH rota- 
tional transitions toward five galaxies, thus considerably 
expanding the existing database. 

2. Observations 

All observations were carried out with the IRAM 30- 
m antenna. The measurements were made with SIS re- 
ceivers tuned to a single sideband with image sideband 
rejections of typically 7dB. An error in the determina- 
tion of the sideband rejection of ±1 dB would result in 
an additional calibration uncertainty of ±5%. The tem- 
perature scale used throughout this paper is main-beam 
brightness temperature (Tmb)- System temperatures were 
between 540 K and 580 K Tmb in the 3 mm wavelength 
range, between 600 K and 900 Tmb at 2 mm wavelength 
and between 1200 K and 1800 K T M b at 1.3 mm. All tran- 
sitions were observed with either a 512 x 1 MHz niterbank 
or a 864 channel acousto-optical spectrometer (AOS) with 
a total bandwidth of 500 MHz. 

The observed transitions, including their designations, 
frequencies and energies above ground (Mangum & Woot- 
ten 1993, Anderson et al. 1990), and the beamwidth of the 
telescope, are summarized in Table 0. 

For broad lines encountered in external galaxies, it is 
not only the system noise that determines the accuracy 
of line intensies and ratios, but also the flatness of the 
spectrometer bandbass, which is affected by the stability 
of the weather and the electronic components. We used 
a wobbling secondary mirror, switching with a frequency 
of 0.5 Hz between the source and two reference positions 



Table 1. Summary of the observed transitions. For H2CO, o 
and p denote ortho and para transitions. 



Line 


Frequency 









GHz 


K 


II 


H 2 CO 








2l2 — In 


140.840 


17 


16 


2o2 — loi p 


145.603 


10 


16 


2n — lio 


150.498 


23 


16 


3l3 — 2i2 


211.211 


32 


12 


3o3 — 2o2 p 


218.222 


21 


12 


3i2 — 2n 


225.698 


33 


12 


CH30H 








5-i - 4 E 


84.521 


32 


27 


2 k - l k 


96.741 


5-20 


27 


3i - 2iA+ 


143.866 


28 


16 


3 k - 2^ 


145.103 


12-51 


16 


7i - 80E 


220.07 


88 


12 


8-1 — 7oE 


229.75 


81 


12 


3-2 -4_iE 


230.027 


32 


12 



a) Energy of the upper level; for E-type methanol, the values 
refer to the li-level, 6.9 K above ground. 

b) Lines contributing to the 2k — lk group: 2_i — 1~\E (5K), 
2 - 1 A+ (6K), 2 - loE (12K), 2i - UE (20 K) 

c) Lines contributing to the 3k — 2k group: 3o — 2qE (19 K), 
3-i -2-iE (12K), 3o -2 A+ {UK), 3 2 - 2 2 A~ (51 K), 

3 2 - 2 2 E (28K), 3_2 - 2- 2 E (32K), 3i - 2 X E (27K), 
3 2 - 2 2 ^1 + (51 K) 



placed symmetrically at 4' offset in azimuth. The back- 
ends are of high stability. Thus, the baselines we obtained 
are very flat. Therefore, only baselines of order zero or 
one were subtracted from the spectra, and the data were 
smoothed to a resolution of ~ lOkms" 1 . 

From continuum cross scans on nearby sources, made 
every ~ 3 hours, we estimate the pointing accuracy to 
be of order 5" or better. To check the calibration of the 
telescope, we usually observed line emission from strong 
Galactic sources like W3(OH) in each period to achieve in- 
ternal consistency. We expect the calibration to be correct 
on a level of ~ ±15%. 



3. Results 

3.1. Formaldehyde 

The parameters of Gaussian fits to the detected lines are 
given in Table|^. Spectra are displayed in Fig. |l|. Since not 
all lineshapes are exactly Gaussian, we have also deter- 
mined the integrated line intensities by adding the chan- 
nels with emission, calculating the errors from the rms 
noise per channel. The results are identical, with the er- 
rors of the Gaussian typically being slightly larger than 
the errors obtained using the rms noise. 
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Fig. 1. H2CO spectra for those positions and transitions where we detected, at least tentatively, a line. All spectra have 
been smoothed to a resolution of lOkms -1 . The central positions (1950.0) and distances of the galaxies are: NGC 253: 
a = 00 h 45 m 06.0 s , 8 = -25°33'36", 2.5 Mpc, Maffei2: a = 02 h 38 m 08.5 s , 8 = 59°23'24", 5Mpc, IC342: a = 03 h 41 m 57.5 s , 
8 = 67°56'25", 1.8 Mpc, M82: a = 09 h 51 m 43.0 s , 8 = 69°55'00", 3.3 Mpc, NGC 6946: a = 20 h 33 m 48.0 soc , 8 = 59°59'00', 11 Mpc. 



We detected mm wavelength H 2 CO transitions toward 
all galaxies that were observed, although the detection of 
the 2o2 — loi line in NGC 6946 is tentative. 

Toward NGC 253, four of the five transitions observed 
could be detected; in M82, all five lines were found. The 
three positions observed in M 82 are detected in the 2i2 — 
In line, while the 2n — lio line is not seen toward the 
(10", 10") offset. Toward this position, the J = 2 - 1 twin 
lines thus have different intensities. 



In IC342, four out of six transitions were detected. It 
is interesting to note that the 3i3 — 2i2 transition was de- 
tected both toward the central and (tentatively) toward 
an offset position, while its twin transition, the 3i2 — 2n 
line, only observed toward the offset position, was not de- 
tected at an rms level of <~ 12 mK. As in M82, the inten- 
sities of the J = 2 — 1 twin lines differ. Toward the center, 
the 2n — lio line is considerably stronger than its coun- 
terpart, the 2i2 — In transition, which is not detected. 
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Table 2. Parameters of H2CO lines observed toward nearby 
spiral galaxies. The offsets from the central position are given 
in arcseconds. Errors (la) are given in parenthesis. 



Source 


Trans. 


/ T M Bdv 


VLSR 


Avi/2 Tmb 








Kkms" 1 


kms 1 mK 


NGC253 












(0,0) 


2(12 


,a),b) 
A 01 


~ 10 


210 


170 ~45 




3o3 


— 2o2 


< 2.8 C) 




- < 17 d) 




2n 


-ho 


6.1(0.7) 


206(7) 


204(25) 28(10) 




3l3 


-2i 2 


5.2(0.7) 205(10) 


151(20) 32(11) 


(-5-2) 


2-12 


- Ill 


8.0(0.9) 


267(6) 


104(13) 72(13) 


(-10-10) 


2u 


-lio 


5.6(0.8) 


284(7) 


126(23) 42(11) 




3l3 


-2l2 


5.5(0.7) 


281(8) 


136(26) 38(11) 


Maffei2 












(0,0) 


3o3 


— 2q2 


< 4.9 


- 


<33 




2l2 


- Ill 


1.7(0.4) 
1.6(0.3) 


-87(6) 
10(6) 


74(18) 22(6) 
69(17) 21(6) 




2n 


-lio 


0.9(0.2) 
0.9(0.2) 


-75(4) 
12(5) 


47(15) 19(4) 
45(11) 18(4) 




3l3 


- 2i2 


< 1.8 




<12 




3l2 


- 2 n 


< 1.6 




<H 


1G342 












(0,0) 


2o2 


— 1 ' 


1.0(0.1) 


27(2) 


47(5) 20(3) 




2l2 


- Ill 


< 0.9 




<16 




2n 


- lio 


1.2(0.2) 


29(4) 


47(11) 23(5) 




3l3 


- 2i2 


1.0(0.2) 


32(6) 


48(12) 19(7) 


Cn <\ 

(U,t>J 


3o3 


— 2o2 


< 3.3 




— <-4U 






^02 


— J-01 


0.6(0.2) 


38(5) 


£oyoj zo^iuj 




2io 


— In 

J- 11 


1.5(0.2) 


44(2) 


28(M 4-9(7) 




^11 


J- 10 


0.3(0.2) c) 


50(2) 


1 1(3") 9Q(7) 

ill 1 A>y\l J 




*J13 


^12 


0.4(0. 2) c) 


27(3) 






Ol2 


^11 


< 0.9 




^1Z 


M82 












(U,U) 


2(12 


1 a ) 


< 3.0 




/on 
- <30 




^12 


— til 


5.3(0.9) 179(17) 


183(37) 27(9) 




2u 


-lio 


2.5(0.6) 191(18) 


147(44) 16(9) 




3l3 


-2l2 


1.3(0.4) e) 


219(16) 


66(25) 16(10) 


(-10,0) 


3o3 




6.6(0.5) 


180(6) 


165(12) 38(8) 


(-10,-10) 


2o2 


— J-oi 


~ 1.8 


120 


110 ~ 15 




2l2 


- Ill 


2.2(0.5) 128(18) 


133(31) 16(10) 




2xi 


-lio 


2.9(0.6) 127(10) 


114(34) 23(11) 




3l3 


-2 12 


3.0(0.5) 


117(8) 


71(14) 40(11) 


(10,10) 


2l2 


- Ill 


5.3(0.6) 


304(5) 


79(10) 63(13) 




2xi 


-lio 


< 0.8 




- <10 



NGC6946 
(10,-10) 2 02 



lox 



l.l(0.5) c) 96(18) 85(22) 13(11) 



a) In NGC253 and M82, the 2 , 2 - lo,i-line is blended with 
the 16-15 transition of HC 3 N, offset by 42 MHz (86.5 km s" 1 ). 
An approximate fit for NGC 253 was obtained by fixing the 
central velocity and width of the HC3N line to the values 
known from the 15-14 and 17-16 HC3N transitions 
(Mauersberger et al. 1990). For M82, the H 2 CO line 
parameters were fixed to the mean value derived from 
unblended H2CO lines in the same position, assuming the 
nominal offset for HC3N. Both estimates are indicated by 
italics. 

b) Data from Mauersberger et al. 1995 

c) Limits for the integrated intensity are 3<r limits, based on 
the rms per channel extrapolated to a linewidth fixed to the 
value determined from detected H2CO transitions toward the 
same position. 

d) All limits and errors to Tmb are rms in a 10 kms -1 wide 
channel 

e) tentative 

f) Data from Baan et al. 1990 



Toward a northern offset position however, the 2i2 — In 
line is stronger than the (tentatively detected) 2n — lio 
line. 

Toward Maffei2, two out of four lines were found. 
The profiles show a double peaked structure. This dou- 
ble feature is also visible in the HNC(l-O) transition 
(Hiittemeister et al. 1995), but is not seen in other molec- 
ular species. 

3.2. Methanol 

The results of Gaussian fits for all observations of single 
positions are given in Table ^ spectra are displayed in Fig. 
pi We also obtained a map of NGC 253 in the Jk = 3k — 2k 
transitions at 145 GHz with a 10" spacing (Fig.|3|). 
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Fig. 2. CH3OH transitions detected toward NGC 253, Maffei 2, 
IC342 and NGC 6964. 



As is the case for the J = 2k— lk transitions at 96 GHz, 
the J = 3k — 2k 'line' is actually a blend of several indi- 
vidual transitions (see Table |l|) with frequencies ranging 
from 145.09375 GHz to 145.13346 GHz, corresponding to 
a velocity shift of 82 kms - . We adopted a rest frequency 
of 145103. 23, GHz, corresponding to the 3 — 2 A+ tran- 
sition. The 3k — 2k transitions were detected with relative 
ease in all sources we observed, with the notable excep- 
tion of M82. In Maffei 2, the 2k — lk line was also easily 
detected. The 5_i — AqE transition in NGC 253 is clearly 
detected. 
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Fig. 3. Map of CH3OH (3k — 2k) toward NGC 253. The upper left panel shows the spectra with a ulsr range of 50-450 kms 
and Tmb ranging from -70 to 200 mK. The total integrated intensity map is given in the upper left panel; solid contours start 
at 5.25Kkms _1 (3a) in steps of 1.75Kkms _1 (la). The dashed contour denotes the 2a level at 3.5Kkms _1 . In the lower 
panel, we present channel maps of 50 kms -1 width. The first solid contour is at 1.8Kkms _1 (3a), the steps are 0.9 K kms -1 ; 
the dashed contour is 1.2Kkms _1 (2a). 



4. Discussion 

4-1. Physical conditions 
4.1.1. H 2 CO 

Formaldehyde is a slightly asymmetric top molecule that 
exists in ortho- (K a = 1,3,5...) and para- (K a = 
0, 2, 4. . . ) forms that are not connected by allowed transi- 
tions. Thus, they can be considered independent molecu- 



lar species. The asymmetry splits each rotational level in 
the K a > 0-ladders into a so-called if -doublet (see, e.g., 
Henkel (1980) for a level diagram and more details). 

Since relative populations of different K a ladders are 
determined by collisions, line intensity ratios of transitions 
from different K a ladders of a given species are excellent 
tracers of the gas kinetic temperature (Mangum & Woot- 
ten 1993). The level populations within a K a ladder are 
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Table 3. Parameters of CH3OH lines observed toward nearby 
spiral galaxies. The offsets from the central position are given 
in arcseconds. 



Source 


Trans. 


J TmbG^ 


ULSR 


Afl/2 


1MB 








Kkms" 1 


kms 1 


mK 


NGC253 














(0,0) 


fi 1 

u— 1 


— 4n F 


2.2(0.4) 


176(6) 


70(12) 


29(9) 




9, _ 


1 a ) 


o.o^u.u 1 


1 f\C\( K\ 

iyu(o ) 


7A/ 1 o\ 

f\J{Lz ) 










14.9(0.7) 


OKJiyyZ ) 


1 f\Af f\\ 
1U4^D ) 


136(15) 




a, 

01 — 


9, 4 + 


< 9 7 b > 






< 9fV 1 Si 




Q 

Ok - 


2k 


21.4(1.8) 


91 1(8.} 




1 02 ('251 




8 1 

O— 1 


— 7nF 


< 9 n 










3-2 


- 4_i_B 


< 5.4 






<19 


Maffei2 














(0,0) 


0—1 


— 4n K 


< 1.6 






<12 




2k - 


lk 


2.8(0.7) 


1 Of 1 f\\ 


77 f 


34(121 




01 




< 0.6 






<6 




3k - 


2v 


4 4f0 71 


1U ) 


1 OOfOf\~\ 


34(111 






7i R 


< 3.9 






<39 




8-1 


- 7o-E 


< 2.7 






<27 


IC342 














(0,0) 


-j— 1 


— 4n F 


< 1.1 






<16 




3 

2 


— 4 1 F 


< 2.1 






<31 


(0,-5) 


3 k - 


2 k 


2.0(0.4) 


33(6) 


48(9) 


40(13) 


(0 5) 


^k 


1 a) 
L k 


4.6(0.3) 


AdfO\ 


Ooyo ) 


68(15) 




Ok 


9, 
^k 


3.9(0.6) 


oy(o ) 


KCtfQ\ 


73(14) 


(0 15) 


0— 1 


— 4n F 


< 0.7 






<11 




Ok 


9, 
^-k 


2.0(0.5) 


oo^D ) 


ARfl 9^ 


41fl3l 






7i F 
* 1 ^ 


< 2.6 






<39 




3_ 2 


- 4_ 1 i? 


< 1.5 






<22 


M82 














(0,0) 


5_i 


-4 E 


< 1.1 






<10 




3 k - 


2 k 


< 2.3 






<15 




80 - 


7iE 


< 2.4 






<18 




3-2 




< 1.6 






<H 


(-10,0) 


^k 


1 a) 


< 3.3 






<34 


f— 10 -10) 


5-i 


— 4o E 


< 1.1 






<10 




3 k - 


2 k 


2^ l.o 






^11 
^11 




80- 


7i-B 


< 1-2 






<H 




3-2 


- 4_i£ 


< 1.8 






<16 


(10,10) 


3 k - 


2 k 


< 1-7 






<15 




80 - 




< 2.5 






<25 


NGC 6946 














(4,0) 


2 k - 


1 a) 
L k 


1.7(7) 


55(7) 


100(7) 


17(7) 


(10,-10) 


5-i 


-4 E 


< 0.7 






< 6 




3i - 


2 X A+ 


< 0.8 






<7 




3 k - 


2 k 


2.5(0.3) 


54(7) 


105(15) 


23(5) 




80 - 


7iE 


< 1.5 






<13 




3-2 


- 4-rE 


< 1-4 






<12 















a) Results taken from Henkel et al. 1987. 

b) Limits to the integrated intensity and Tmb were derived as 
for H 2 CO (Table §) 



determined by an equilibrium between collisional excita- 
tion and spontaneous and collisional deexcitation, making 
the relative level populations within a K a ladder sensi- 
tive to the H2 density and, to a lesser degree, to kinetic 
temperatures. 

We have carried out large velocity gradient (LVG) sta- 
tistical equilibrium computations of the population of the 
K a = 1 ladder of ortho formaldehyde (o-H 2 CO), for which 
we have the most complete data. The model is described 
in detail in Henkel et al. (1980). We included 16 levels of 
that rotational ladder fixing the kinetic temperatures to 
a value of 80 K. Setting Tkm to 50 K changes the line in- 
tensity ratios only by a few percent. As a result of these 
calculations, line ratios for the 140, 150 and 211 GHz tran- 
sitions of 0-H2CO are related to the H2 density and the 
H2CO column density per unit linewidth (see Fig.||). 

It is evident from Fig.^ that the 1.3 mm and 2 mm 
lines have similar intensities if densities exceed 10 5 cm -3 
(this also holds for para formaldehyde). In order to de- 
rive column densities of 0-H2CO from the integrated in- 
tensities of one or several lines, an assumption must be 
made about the H2 density and, to a lesser degree, the 
kinetic temperature, both of which determine the rota- 
tional temperature, or more generally, the partition func- 
tion. An integrated intensity of 1 Kkms -1 of the 140 GHz 
line corresponds e.g. to Af(o — H2CO) ~ 510 13 cm~ 2 if 
one assumes a low H2 density of 10 4 cm~ 3 ; if, however, 
n(H 2 ) = 10 5 — 10 7 cm~ 3 is assumed, then the o-H 2 CO 
column densities for the same line intensity are an order 
of magnitude lower. In Table[|, we estimate iV(o — H2CO) 
from the intensities of the 140 GHz lines, and the den- 
sities from the measured line intensity ratios. The most 
reliable densities from mm-wave line intensity ratios are 
estimated for regions where n(H2) > 10 5 cm -3 . For lower 
densities, mm-wave lines of H2CO become optically thick 
if 7V(o-H 2 CO)/Aw £ 10 12 - 5 cm- 2 /kms -1 . In this case, 
the line intensity ratios not only depend on n(H2), but also 
strongly on iV(o — H2CO)/Aw. Upper limits on n(H2) es- 
timated in Tabled are based on the optically thin limit. 

In order to determine relative abundances, the H2 col- 
umn density has to be known. In Table ||, we use values de- 
duced from CO measurements, and the standard conver- 
sion formula of Strong et al. (1988), i.e. 77AT(H 2 )/J(CO) = 
2.310 20 cm~ 2 (Kkms -1 ) -1 . However, there is growing ev- 
idence that this conventional conversion factor cannot 
be applied to Galactic bulge regions (see Dahmen et 
al. 1996a, b for our Galactic center, Mauersberger et al. 
1996a,b for NGC 253 and NGC 4945, Shier et al. 1994 for 
IR- luminous galaxies) . The correction factor 7/ might be as 
large as 10, and is likely to be different for different galax- 
ies. It will, however, always decrease H2 column densities 
and increase the relative abundance of the molecule. 

NGC 253. Toward this galaxy, we have the largest data 
base. In order to compare line intensities which have been 
obtained at different wavelengths one has to take into ac- 
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count the different beam sizes and the spatial structure of 
the source. Since we did not map the distribution of H2CO 
toward any of our sources, we assume that the distribu- 
tion of this molecule toward NGC 253 is similar to that 
of the J = 2 — 1 line of 12 CO, which has been mapped 
with 12" resolution and convolved to different resolutions 
by Mauersberger et al. (1996a). From an interpolation of 
the data toward the central position of NGC 253 in their 
Table 1, intensities in a 16" beam are a factor 0.85 weaker 
than data at 12" resolution. In the following analysis we 
assume the same correction factor for the 2 mm data of 
M 82 and IC 342 where the size of the molecular emission 
is similar to that of NGC 253. 

We can use the intensity ratios of lines within the 
K a = 1 ladder of ortho formaldehyde to estimate the 
H2 density of the H2CO emitting gas. In practice, a dif- 
ficulty in assessing the line ratios arises if the line shapes 
and/or the velocity coverage differs. One might then in- 
correctly compare different gas components. The 3i3 — 2 i2 
(211 GHz) and the 2 12 - In (140GHz) lines were ob- 
served toward slightly different positions, adding to the 
uncertainty. If only the clear overlap region in the spec- 
tra is considered, a line intensity ratio of 0.42 (±0.07) re- 
sults. From Fig.H such a low ratio is compatible with 
n(H2) ~ 10 4 cm" , and even lower densities if H2CO lines 
are saturated. We conclude that n(H 2 ) < 10 5 cm -3 for 
the bulk of the H2CO emitting gas. From our limit to 
the 3o3 — 2o2 line and the detection of the 2o2 — lot line 
of para H2CO, the corresponding ratio is < 0.2 which is 
compatible with n(H2) < 10 4 cm~ 3 . 

The ratio of the 2n — lio (ortho, 150 GHz) and the 
2o2 — loi (para, 146 GHz) transitions depends on the or- 
tho/para (o/p) ratio. For the high temperature limit of the 
o/p ratio of 3, one would expect a line intensity ratio of 
- 1.4 (for T kin = 80 K) and n(H 2 ) = 10 4 c m - 3 . The value 
we actually observe is 0.6, or even lower, if one argues 
that not the entire velocity range for the wide 2n — lio 
line should be considered for the ratio. We therefore sug- 
gest that the o/p ratio of H2CO actually has a value close 
to 1, which is lower than the value derived by Aalto et al. 
(1997) for the 1 mm transitions. Since the o/p ratio bears 
important information on the formation of molecular gas, 
and our result, involving a blended line, is uncertain, fur- 
ther observations are very desirable. 

M 82. We have measured ortho-H2CO lines toward two 
positions in M82, namely toward the nucleus and to- 
ward the so called SW molecular hotspot at an offset 
of (Aa,AS) = (-10", -10") (Mauersberger & Henkel 
1991). The ratio (for the same beam size) of the 3i3 — 2i2 
(211 GHz) and the 2 12 - In (140 GHz) lines is < 0.2 
toward the nucleus, regarding the tentatively detected 
3i3 — 2 12 as contributing an upper limit only. Toward the 
SW hotspot, however, it is as high as 1.2(±0.4), or even 
~ 1.5, if we restrict the intensity from the 2 12 — In line 
to the width of the slightly narrower 3i3 — 2i2 transition. 




4 5 6 7 

og n(H 2 ) 

Fig. 4. Predicted line intensity ratios of ortho-formaldehyde 
lines as a function of n(H2) (in cm -3 ) and jV(ortho-H2CO) /An 
(in cm" 2 /(kms _1 ). Solid contours: T B (211 GHz)/T B 140 GHz), 
thin contours: T B (150 GHz)/Te(140 GHz). The assumed T kin 
is 80 K. 

Thus, toward the nucleus, n(H2) < 10 4 cm~ 3 ; while den- 
sities are as high as n(H 2 ) ~ 10 6 cm~ 3 toward the SW 
hotspot. This estimate is relatively robust to whether or 
not the H2CO lines are saturated. 

Such density variations were already predicted from an 
analysis of centimeter wave transitions of ortho formalde- 
hyde observed with lower angular resolution by Baan et 
al. (1990). While the 6 cm line shows a broad absorption 
(Graham et al. 1978) over the entire velocity range, the 
2 cm line is seen in emission over a velocity range that hints 
toward an origin which is confined toward the southwest- 
ern part of the nuclear region. The 2 cm H2CO line is usu- 
ally observed in absorption, even against the 2.7 K cosmic 
background. 2 cm emission requires H2 densities exceeding 
10 5 ' 5 cm -3 . This is in very good agreement with our re- 
sults from millimetric lines toward the SW hotspot. Also 
the central concentration of N2H" 1 ", a molecular ion that is 
destroyed in a high density medium, indicates that the gas 
toward the central region has a much lower density than 
toward the SW hotspot (Mauersberger & Henkel 1991). 

From the observed line intensity ratio of the 2i2 — In 
and 2o2 — 2oi lines of ~ 1.2, the o/p ratio toward the SW 
hotspot is (assuming n(H2) = 10 6 cm~ 3 ) ~ 2, i.e. close 
to the high temperature limit. Note, however, that, as for 
NGC 253, this result may be affected by the uncertainty 
of the fit to the blended 2o2 — 2oi transition. 

IC342. This nearby (~ 2Mpc; McCall 1989, Karachent- 
sev & Tikhonov 1993) and nearly face-on spiral galaxy is 
very similar to the Milky Way galaxy with respect to the 
IR luminosity and the gas mass in its central region. 

We have observed formaldehyde toward the nucleus 
(0",0"), where the velocity coverage and the line shapes 
of all transition are in excellent agreement, and toward one 
offset position, (0", 15"), where the (tentative) 3i3 — 2i2 
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transition seems to appear at a different velocity. The 
beam toward the (0", 0") position also contains the molec- 
ular clump B (Downes et al. 1992), which is associated 
with free-free emission equivalent to 300 05 stars or 30 
times more than the emission from the Galactic center 
star forming region SgrB2. 

3i 3 - 2 12 /2 U - li (211 GHz/150 GHz) line ra- 
tios indicate high densities of 10 62 cm" 3 toward the 
(0",0") position. The limit to the 2 U - lio/2i 2 - In 
(140 GHz/150 GHz) line ratio toward the nucleus of > 1 is 
also suggesting high densities, while the limit toward the 
(0", 15") offset (< 0.4, considering the tentative detection 
to be an upper limit, and regarding a common velocity 
interval for both lines), indicates a far lower density of 
< 10 4 cm" 3 . As in M82, our two beams pick up contri- 
butions from two molecular phases, namely from a mod- 
erately dense interclump gas and from the clumps them- 
selves. The presence of several molecular gas components 
has also been inferred by Downes et al. (1992) from the 
comparison of single dish and interferometric observations 
in CO. 

Maffei2. The ratios of the 150 and 140 GHz line, which 
agree very well in velocity coverage and line shape, indi- 
cate typical densities of the order of 10 4 cm" 3 . 



group of transitions to determine the CH3OH column den- 
sities given in Table ^. 
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Fig. 5. Predicted intensity ratios of methanol lines for a 
raster of values for n(H 2 ) (in cm" 3 ) and iV(CH 3 OH))/Av 
(in cm" 2 /(kms~ 1 ). Solid contours: T B (145 GHz)/T B (96 GHz), 
dashed contours: T B (145 GHz)/T B (84 GHz). The assumed T kin 
is 80 K. 



4.1.2. CH3OH 



Methanol is an asymmetric top molecule capable of hin- 
dered internal rotation. Transitions from A-type (internal 
rotation) to E-type (no internal rotation) levels are strictly 
forbidden; the two types behave as separate molecular 
species in the ISM. In addition, a and 6-type transitions 
have to be distinguished according the component of the 
dipole moment that is is parallel to the axis of rotation. 
All transitions contributing to the 3k — 2k and 2k — lk lines 
are a-type (Afc = 0; /i a = 0.896 Debye), while the 5_i — 4 
transition is &-type (Afc = ljfib = 1.412 Debye) (Lees et 
al. 1973, Sastry et al. 1981; for a detailed discussion of the 
CH3OH molecule, see e.g. Menten 1987). 

An analysis of the 3k — 2k group of transitions is diffi- 
cult since the lines are blended. We have, however, carried 
out LVG calculations analogous to those for H2CO, using 
a program provided by CM. Walmsley (see, e.g. Walmsley 
et al. 1988, Bachiller et al. 1995). We have modelled the 
intensity ratios of the sum of all transitions contributing 
to the 3k — 2k group to the sum of the 2k — lk transitions 
and the ratio of the 3k — 2k group to the 5_i — line 
as a function of H2 density and CH3OH column density 
per unit linewidth (Fig. ||), for a kinetic temperature of 
80K. Clearly, for 7V(CH 3 OH)/Au ^ 10 14 cm^/kms- 1 , 
the line ratios are independent of column density. Thus, 
the CH3OH lines are very likely to be optically thin in 
all our sources. We have used the H2 densities determined 
from Fig. ^ and the integrated intensities of the 3k — 2k 



Table 4. Beam averaged H2CO and CH3OH column densities 
and relative abundances. 



Source r/N 



N 



N 



[H2CO] 
[H 2 ]ij 



[CH3OH] 

[H2]i) 





H 2 a) 


H 2 CO 


CH3OH 








10 22 


10 13 

cm -2 


10 13 


10- 9 


10" 9 


NGC253 












(0,0) 


37 


40 b) 


50 b > 


1.1 


1.4 


M82 












(0,0) 


16 


~ 25 c) 


<2 e) 


1.6 


< 0.1 


(-10,-10) 




0.9 d) 


< l-2 e) 






IC342 












(0,0) 


4.9 


0.2 d) 




0.04 




(0,5) 


4.9 




u /) 




2.2 


(0,-5) 






6 /) 






(0,15) 




> 8 C) 


6 /) 






Maffei 2 


1.4 


15 6 > 


4&) 


11 


2.9 


NGC 6946 


2.3 




3 i,) 




1.3 



a) based on /(CO) (see Hiittemeister et al. 1995 and 
references therein) and a "conventional" conversion factor. 

b) assumed n(H2): ~ 10 4 cm" 3 

c) assumed n(H2): < 10 4 cm -3 

d) assumed n(H2): ~ 10 6 cm -3 

e) from Eq. 1; with T rot =10 K 
/) assumed n(H 2 ): k, 10 3 cm" 3 
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NGC253. The map we obtained in the 3k — 2k transi- 
tion toward NGC 253 (Fig. |J) shows the same general 
shape as maps in other high density tracing molecules like 
HNC (Hiittemeister et al. 1995) and CS (Mauersberger & 
Henkel 1989). The emission extends from the north-east 
to the south-west, with a (not deconvolved) full source size 
of ~ 55" x 25", which is similar to the extent of intense 
12 CO (2-1) emission (Mauersberger et al. 1996a). 

Therefore, it seems reasonable to correct the 3k — 
2k/2k — Ik and the 3k — 2k/5_i — A E line ratios for the 
different beamsizes in the same way as for H2CO (Section 
4.1.1), based on the Mauersberger et al. (1996a) data. The 
line intensities of the 3 mm lines then have to be multi- 
plied by ~ 1.4 to make them directly comparable to the 
2 mm lines. The linewidths of the two groups of transitions 
agree well. The resulting 3k — 2k/2k — lk ratio (0.8±0.15) 
implies a H2 density of ~ 10 35 cm -3 . The 5_i — 4o line 
is much narrower than the 3k — 2k group. Since we are 
comparing a group of transitions spaced by > 80kms _1 
to an unblended line, this is to be expected. However, it 
might not account for the entire difference. The ratio of 
the integrated intensities (6.8±2.4) thus is an upper limit, 
yielding a lower limit to the density of 10 45 cm -3 . If we 
consider only the velocity interval for the 3k — 2k group 
where 5_i — 4o emission is present, this lower limit to the 
intensity ratio of 3.1±0.7 corresponds to a an upper limit 
to the H2 density of > 10 5 cm -3 (see Fig. ||). In any case, 
the values from the 3k — 2k/2k — lk and the 3k — 2k/5_i — 4 
line ratios embrace the H2 density estimated from H2CO. 
The critical density of the 5_i — AqE transition is smaller 
than the mean of the 3k — 2k group by a factor of ~ 5, 
while the difference is much less pronounced between the 
3k — 2k and the 2k — lk groups. This might explain the 
difference in the density estimates resulting from the two 
line intensity ratios. 

At an H2 density close to 10 4 cm -3 , both the 2k — lk 
and the 3k — 2k groups of lines are subthermally excited 
(T C x< 7ki n ), with excitation temperatures for the individ- 
ual transitions ranging from 10 K to 20 K. We thus obtain 
a value for iV(CH30H) which is lower by a factor of ~ 6 
than what is determined under the assumption that the 
density is sufficiently high to thermalize the CH3OH emis- 
sion (Henkel et al. 1987). 



M 82. The clear non-detection of any methanol line in 
M82 is one of our most striking results. Henkel et al. 
(1987) searched for the 2k — lk transition in this galaxy 
and established an upper limit of 30 mK rms Tmb for 
this line. We obtained upper limits of 10-15 mK rms for 
the 3k — 2k transition toward the central position and 
the north-eastern and south-western peak positions of the 
molecular ring. 



Since no line ratios could be determined, we have used 
the LTE approximation for optically thin lines to estimate 
the total CH3OH column density: 

1 67 1 14 r 

7V(CH 3 OH) = 1.28 T™ e £ '/ fcT - X '°J ™ J T MB dv (1) 

(E\: Energy of lower level, v. line frequency in GHz, fi: 
dipole moment in Debye, S: linestrength; Menten et al. 
1988). 

Using Eq.[l], we have fitted all components of the 3k— 2k 
group to the limit to the integrated intensity. For subther- 
mally excited lines (T rot = 10 K), we find column densi- 
ties ranging from smaller than 1 10 13 cm -2 (SW hotspot) 
to less than 210 13 cm -2 (center and NE hotspot), corre- 
sponding to an abundance of [CH 3 OH]/77[H 2 ]< 1.210" 10 
(see Table @). If we assume that the H2 densities are the 
same as the ones found for H2CO and apply our LVG 
model, the only change is a slight increase in the limit for 
iV(CH 3 OH) to 2 10 13 cm" 2 in the dense SW hotspot (see 
Table §. 

IC 342. Within the 3k — 2k blend, two groups of lines 
are separated by ~ 45kms . The group of 5 lines at 
higher frequency has generally higher energies than the 
lower frequency group. Because the lines in IC 342 are nar- 
row, these two groups are resolved, and the higher energy 
transitions are below the detection limit. Thus, LTE fits 
to the 3k — 2k group already greatly constrain the possible 
excitation conditions: T ro t has to be 10 K or lower. Since 
the kinetic temperature of the gas in IC 342 is likely to 
be at least 50 K (Ho et al. 1990), CH 3 OH has to be very 
subthermally excited, with T ex lower than for NGC 253. 

The LVG calculations confirm this: For a (beam size 
corrected) 3k — 2k/2k — lk intensity ratio of 0.6(±0.1), ob- 
served toward the (0",5") position, we find an H2 density 
of only - 10 3 cm- 3 (T kln = 80 K) to - 10 3 4 cm- 3 (T kin = 
50 K) and excitation temperatures of ~ 5 — 10 K. This im- 
plies that CH3OH is sensitive to interclump gas of only 
low to moderate density. 

The beam-averaged column density we derive is 
1.1 10 14 cm~ 2 for a position slightly north of the center of 
IC 342, and half that for offsets to the south and further to 
the north (assuming the same low density), corresponding 
to an abundance of [CH 3 OH]/7?[H 2 ] ~ 2 10" 9 (Table @). 
We note that the column densities given in Henkel et al. 
(1988), based on T rot = 50 K and derived from the 2k — lk 
line, are far too high. A recalculation using T rot = 8 K 
and Eq. [I] (with ^ Si = 7 and Ei = 6 K) yields a total 
column density that agrees to within 10% with our result 
for the 3k — 2k line, showing that the LTE and LVG ap- 
proximations converge in this case. The limit given by the 
non-detection of the 5_i — 4q transition is also compatible 
with a total CH3OH column density of ~ 10 14 cm -2 . 

Maffei2 and NGC 6946. The 3 k - 2 k /2 k - l k line ratios, 
corrected for beam size, imply H2 densities of ~ 10 4 cm~ 3 . 
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It thus seems likely that the excitation of CH3OH in both 
galaxies is similar to what is found in NGC 253. 

4-2. Chemical variations 

Molecular abundances. For most of the galaxies observed, 
the column densities of o-H 2 CO and CH 3 OH are compara- 
ble. This might also be the case for the M82 SW hotspot, 
where we have only a limit for the CH3OH column den- 
sity. The notable exception is the nucleus of M 82 where 
the column density of CH3OH is at least an order of mag- 
nitude lower than that of H2CO. It is quite clear that this 
is not just mimicked by excitation effects. Such a discrep- 
ancy between M 82 and other starburst galaxies, especially 
NGC 253, has already been noticed before (e.g. Mauers- 
berger & Henkel 1993). 

Evaporation from grain surfaces is in many cases a 
crucial process to maintain a high gas phase abundance 
of complex molecular species. Models predict that these 
molecules, such as CH 3 OH and also H 2 CO, are chemi- 
cally converted into simpler species on a timescale of or- 
der 10 5 years (Helmich 1996). Methanol evaporates at a 
temperature around 70 K (Turner 1989, Nakagawa 1990) 
and is known to be a tracer of hot, dense gas (Menten 
et al. 1988), since its abundance is observed to increase 
dramatically in the vicinity of young massive stars. 

Formaldehyde, on the other hand, empirically shows 
far smaller abundance variations in Galactic interstellar 
clouds (Mangum and Wootten 1993) although it is 3- 
25 times more abundant in the Orion hot core than in 
other molecular clouds (Mangum et al. 1990). Contrary to 
CH3OH, H2CO is also abundant in cirrus clouds (Turner 
1993) and cool, dense Galactic disk sources. It seems cer- 
tain that processing on dust grains must play a role in 
the chemistry of H2CO, but the replenishment process is 
unclear (Federman & Allen 1991, Turner 1993, Liszt & 
Lucas 1995). 

Takano et al. (1995) notice that all molecules that are 
known to be depleted in M82 (besides CH 3 OH also SiO, 
HNCO, CH 3 CN and SO) form preferentially under high 
temperature conditions. In the presence of a steep grav- 
itational potential and a central bar, we expect shocks 
and turbulence to play a significant role in the heating of 
molecular clouds distant from the actual sites of star for- 
mation (e.g. Hiittemeister et al. 1993, Das & Jog 1995). 
This can explain the high temperatures toward the central 
molecular condensation of NGC 253, which is much more 
compact than that of M82, a smaller, irregular galaxy. 
Tidal heating should therefore operate less efficiently in 
M82 than in larger spirals like NGC 253, NGC 6946 or 
IC 342. Thus, one would expect the bulk of the gas in the 
center of NGC 253 to be warmer than that toward the 
center of M 82 (with the exception of those clouds which 
are heated directly by newly formed stars). 

The abundance of CH3OH in NGC 253 does not seem 
to be exceptionally high when compared to NGC 6946 or 



non-star burst nuclei like Maffei2 and IC342. But we ten- 
tatively find that the gas traced by CH3OH in NGC 253 
is at a similar density as the gas traced by H2CO, while 
the gas traced by CH3OH close to the nucleus of IC 342 
is at a lower density. If temperature is indeed the key 
to the distribution of CH 3 OH, we may conclude that the 
warm gas in IC342 is less dense than in NGC 253. Since 
IC 342 is not a starburst galaxy, this agrees with the fact 
that warm gas in the center of our Galaxy is, in the ab- 
sence of massive star formation, at lower densities than 
cooler gas (Hiittemeister et al. 1993), while the opposite 
is expected for active star forming regions. NGC 253 might 
be lacking the high temperature, low density interclump 
gas component picked up by CH3OH in IC342; the bulk 
of the molecular material in NGC 253 is warm, at least 
moderately dense gas. In this scheme, M82 has a cooler 
interclump component, seen in H2CO and N2H + (Mauers- 
berger & Henkel 1991), but not in CH3OH. 

Extended methanol emission in NGC 253. While the gen- 
eral shape of the central molecular condensation seen in 
CH3OH and other molecules is similar, in a detailed com- 
parison differences are apparent. The maps of total inte- 
grated intensity in 12 CO(2-l), CS and HNC all show only 
one central peak. Only when the red- and blueshifted emis- 
sion is displayed separately, two distinct peaks, separated 
by 10" — 15", become visible, with the redshifted emis- 
sion centered toward the south-west. The CH3OH map 
(Fig. H) shows these two peaks even for the total velocity 
range. The positions of the peaks agree to within better 
than 5" with the 'red' and 'blue' peaks seen in the other 
molecules. From the channel maps, the blueshifted emis- 
sion is strongest at 200 km s _1 and the red-shifted emission 
peaks at 300 km s _ . In addition, there is a third peak, vis- 
ible in the intensity and channel maps, at an offset of ~(- 
25", -10"), at an even higher velocity of 350 - 450 kms -1 . 

Only the interferometric 12 CO(1-0) map obtained by 
Canzian et al. (1988) with a 5" x 9" beam also shows 
two resolved peaks (at the same positions as the CH3OH 
peaks) when the total velocity range is considered. Thus, 
the single dish CH3OH data match the interferometric 
data not only more closely than the single dish 12 CO data, 
but also better than single dish maps of other high den- 
sity tracers, namely CS and HNC. Since the interferometer 
misses extended emission, we can conclude that CH3OH 
traces a more confined gas component than CS or HNC. 
Since the critical densities of CH 3 OH(3 k - 2 k ), CS(2-1) 
and HNC(l-O) are roughly similar at ~ 3 — 610 5 cm~ 3 , 
the differences in spatial distribution seem to be caused 
by chemical fractionation rather than gas density. Since 
production of CH3OH is favored in high temperatures, 
the differences in distribution may reflect differences in 
gas temperatures within the generally at least moderately 
dense gas in the bulge of NGC 253, with CH3OH preferen- 
tially tracing the sites of ongoing massive star formation. 
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5. Conclusions 

We have observed rotational transitions of the high gas 
density tracing molecules H2CO and CH3OH toward 
a sample of nearby gas rich external galaxies, namely 
NGC 253, Maffei 2, IC 342, M 82 and NGC 6946. Our main 
results are: 

1. H 2 CO was detected in all 5 galaxies searched (tenta- 
tively in NGC 6946). 

2. In the prominent starburst galaxy NGC 253, the 
bulk of the gas emitting H2CO lines is at densi- 
ties of ~ 10 4 cnT 3 . In both M82 and IC342, the 
H2CO line ratios distinguish two different compo- 
nents of molecular gas. High densities dominate the 
SW molecular hotspot in M 82 and the central region 
of IC342 (n(H 2 ) ~ 10 6 c m - 3 ), while lower densities 
(< 10 4 cm~ 3 ) are prevalent toward the center of M82 
and the north of the nucleus of IC 342. 

3. We detected CH3OH lines in the same sample of galax- 
ies, with the notable exception of M82. 

4. In all galaxies detected, methanol is subthermally ex- 
cited. In NGC 253, Maffei 2 and NGC 6946, moderate 
gas densities of ~ 10 4 cm~ 3 are likely. Toward a posi- 
tion north of the center of IC 342, the H2 densities seem 
to be at ~ 10 3 cm -3 . This might indicate that warmer 
gas is at lower densities in the non-starburst nucleus of 
IC342, where CH3OH traces an entirely different gas 
component than H 2 CO. 

5. In M82, no methanol could be detected down to an 
rms level of 10-15 mK, corresponding to an abundance 
at least an order of magnitude lower than in the other 
galaxies. This confirms large scale chemical differences 
in starburst centers. The likeliest cause of these are 
global differences in the temperature of the molecular 
material. 

6. A map of NGC 253 in the 3k— 2k line of methanol shows 
several peaks of emission and a greater amount of spa- 
tial structure than than single-dish maps of of other 
molecules, even high-density tracers. Thus, CH3OH 
emission is associated with smaller scale clumps. This 
indicates chemical fractionation not only between dif- 
ferent galaxies but also within the bulge of NGC 253, 
possibly reflecting the temperature structure of the 
gas. 
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